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FOREWORD 


This  report  was  prepared  at  the  Electronic  Sciences 
Laboratory  as  part  of  the  Independent  Research  pro¬ 
gram  of  LMSC .  It  is  the  first  of  a  series  of  contri¬ 
butions  toward  finding  appropriate  diffraction-pattern 
manipulation  to  improve  the  discrimination  probability 
of  aberrant  optical  systems  for  off-axis  objects. 
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ABSTRACT 


This  paper  discusses  a  general  solution  of  the  diffraction  integral  for  a  rectangular 
aperture  considered  in  scanning  procedures  of  diffraction  patterns.  It  is  shown  that 
in  the  presence  of  aberrations,  the  integrand  in  Fresnel-Kirchhoff’s  diffraction  for¬ 
mula  can  be  expanded  into  a  series.  By  this  expansion,  the  two-dimensional  diffrac¬ 
tion  integral  is  converted  into  a  summation,  the  single  terms  of  which  divide  into 
products  of  two  functions  where  only  one-dimensional  integrations  over  products  of 
Legendre  polynomials  are  to  be  performed.  In  these  integrations,  use  is  made  of  the 
orthogonality  property  of  Legendre polynomials. 

Particular  solutions  are  obtained  concerning  the  two-dimensional  diffraction  pattern 
intensity  distribution  in  the  Gaussian  image  plane  for  an  aberrant  optical  system  and 
the  three-dimensional  intensity  distribution  near  the  focus  of  an  aberration-free 
optical  system. 
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Cartesian  coordinates  of  Q  referred  to  system  with  origin  at  Pj 
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Cartesian  coordinates  of  P  referred  to  system  with  origin  at  P* ; 

these  coordinate  systems  have  mutually  parallel  and  equally  oriented 
axes  with  common  Z-axis  along  the  line  [  C  Pj^  ] 


dimensionless  width  and  length  of  rectangular  area  in  the  exit  pupil , 
respectively ,  corresponding  to  aperture 
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=  -X_ 


dimensionless  coordinates  of  point  Q  on  reference  sphere 


u  =  kZ 
V  =  kX 
w  =  kY 


dimensionless  ’’optical  coordinates  ’  of  observation  point  P 


X  ,  Y 
o  o 


X^Y* 


Cartesian  coordinates  of  referred  to  system  with  origin  at  the 
center  of  the  object  plane 

Cartesian  coordinates  of  P*  referred  to  system  with  origin  at  the 
center  of  the  image  plane;  arbitrarily  choosing  the  meridional  plane 
defined  by  optical  axis  and  principal  ray  as  reference  plane ,  one  gets 
X*  =0 

O  1 

Cartesian  coordinates  at  point  P^  where  a  specific  ray  emerging 
from  P^  intersects  the  aperture  plane;  reference  system  has  origin 
at  the  object  plane  center 
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Xj  ,  Yj  ,  Zj  Cartesian  coordinates  of  point  P'j^  in  the  exit  pupil  plane  which  corre¬ 
sponds  to  (with  or  without  diffraction);  reference  system  has  origin 
at  the  image  plane  center 


Xi,Y, 


X’  ,  Y’  ,  Z* 


X*  ,  Y*  ,  Z’ 


Cartesian  coordinates  of  point  in  the  image  plane ,  aberrant  by 
small  amounts  from  X^  and  Yj 

Cartesian  coordinates  of  observation  point  P  in  the  neighborhood  of 
P^  ,  referred  to  system  with  origin  at  the  image  plane  center 

Cartesian  coordinates  of  intersection  points  and  Q  of  rays 
[  P’^  P j  ]  and  [  P 1  >  respectively ,  with  Gaussian  reference  sphere 

centered  at  P^  and  passing  through  the  center  C  of  exit  pupil 

Cartesian  coordinates  of  intersection  points  and  Q  of  rays 
I P^  P^  ]  and  [  P’j^  P  ]  ,  respectively ,  with  actual  wavefront  through 
C^;  the  coordinate  systems  used  here  have  also  mutually  parallel  and 
equally  oriented  axes  with  common  Z^axis  along  the  optical  axis  of  the 
image  forming  instrument  (positive  in  the  sense  of  light -wave  propa¬ 
gation  from  the  object  to  the  image  space) 
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>  Seidel  variables 


refractive  indices  of  object  and  image  spaces »  respectively,  without 
restricting  the  generality  one  may  use  n^  -  n^  ^  1 

distance  between  object  and  entrance  pupil  planes 

distance  between  image  and  exit  pupil  planes;  for  convenience,  «  R 
is  used  where  R  is  the  radius  of  the  Gaussian  reference  sphere 
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o  1  *^1 

^2  >2^2 
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K  =Xili+ypi 


X*  =  -XjX*/R; 


unit  lengths  in  the  entrance  and  exit  pupil  planes ,  respectively , 
is  arbitrarily  chosen  equal  to  unity 

by  Gaussian  optics,  referred  to  Seidel's  coordinates  of 
and  P* 

basic  quantities  needed  for  analytical  description  of  wave 
aberration  in  Seidel  coordinates 


y|;  =  -Vi/^ 

»xAj  =  R^^|/Aj; 


relation  between  Seidel's  coordinates  and  integration 
variables  used 


«  yAi  =  J 

^(2k) 

^(4) 

B 

C 

D 

E 

F 


aberration  function  of  order  (2k  -  1) 
primary  or  third-order  wave  aberration 
coefficient  of  third-order  spherical  aberration 
coefficient  of  third-order  astigmatism 
coefficient  of  third-order  curvature  of  field 
coefficient  of  third-order  distortion 
coefficient  of  third-order  coma 


2.  Symbols  used  in  expansions  and  evaluation 
sin  0,.p  cos  0) 

=  rJJ'(p)  Zemike's  circle  polynomials  in  two  real  variables 

X  =  p  sin  d  ;  y  =  p  cos  0 
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6^  =  1  for  a  =  /3 
a 

=  0  for  a  ^  P 


radial  polynomials ,  n  ,  m  being  arbitrary  indices  restricted 
by  m  I  0  and  n  S  0  integer  ,n2|m|  ,n-|m|  even 

Kronecker  symbol  referred  to  arbitrary  indices  a  ,  p 


Bessel  function  of  order  (s  +  1/2) .  s  being  integer,  with 
arbitary  argument  p 


Pg(cos  e ) 
p  /6' 


K 


I 

n-m  ,  m 


i 


Legendre  polynomial  of  order  s  in  usual  definition 

coefficient  in  defining  the  Legendre  polynomial  of  order  n 
by  a  power  series  ;  6^  is  common  denominator  for  all 
values  ^  is  running  index 

th 

coefficient  in  expanding  the  n*'  power  in  terms  of  Legendre  *s 
polynomials ,  6^  being  again  a  common  denominator  for  all 
values  a 

n ,  V 

constant  in  expressing  the  power  of  the  wave  aberration 
$  as  a  double  series  in  terms  of  power  products 
where  4  i  ^  are  integration  variables 

jp 

index  in  K  like  ( n  -  m )  and  n  ,  as  well  as  running 

n-m  ,  m 

quantity  in  power  series  expansion  of  exponential  function  e^^^ 


,  w )  solution  functions  occurring  in  the  expanded  form  of  the 

disturbance  U( u  ,  v  ,  w )  at  the  observation  point  P( u,  v  ,  w ) ; 
(conversion  of  two-dimensional  diffraction  integral  into  a 
series  containing  only  one-dimensional  integrals) 


p\(x)  = 


new  introduced  polynomials  expressible  as  a  double  series 
of  Legendre  polynomials,  or  even  as  a  single  series  of  the 
form 


= 


r  ,  1 


2r+k-2i 


(X) 
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Gjj(z) 


where  x  is  an  arbitrary  variable,  r,k,i  are  subscripts; 

is  common  denominator  for  all  coefficients  . 
r  r ,  1 

newly  defined  solution  function  of  arbitary  variable  z  , 

obtained  by  applying  orthogonality  property  of  Legendre 

polynomials  to 

1 

/  /^J(x)P^,(x)dx 
-1 


I(u,  V,  w)  =  |U(u  ,v,  w)|^  intensity  distribution  of  diffraction  pattern  at  P(u  ,  v ,  w) 
2 

I^[  ( A/A,)4p^i'^]  central  (maximum)  intensity  of  aberration-free  diffraction 

pattern 
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Section  1 
INTRODUCTION 


Determination  of  the  intensity  distribution  of  a  diffraction  pattern  is  based  on  the  evalu¬ 
ation  of  the  Fresnel -Kirchhoff  diffraction  formula  (Ref.  1).  Although  this  formula  is 
obtained  only  as  an  approximation  from  the  more  general  Helmholtz-Kirchhoff  integral 
theorem,  its  usefulness  and  validity  have  been  proven  in  many  cases  in  the  past. 

Historically,  the  classical  evaluation  technique  of  the  diffraction  integral  yielded  the 
discrimination  between  Fresnel  and  Fraunhofer  diffraction  patterns  which,  of  course, 
have  been  widely  investigated  under  a  variety  of  conditions  for  the  parameters  involved. 
Relatively  soon  after  the  publication  of  the  Fresnel -Kirchhoff  formula  in  1882,  Lommel 
(Ref.  2)  extended  the  evaluation  technique  in  order  to  determine  the  three-dimensional 
intensity  distribution  of  the  diffraction  patterns  in  the  neighborhood  of  the  Gaussian 
focus  of  an  aberration- free  optical  system  due  to  circular  and  rectangular  apertures. 
Later,  attempts  were  made  to  study  the  effect  of  aberrations  on  diffraction  patterns, 
especially  with  regard  to  circular  apertures.  As  pointed  out  by  Wolf  (Ref.  3)  in  a 
historical  review,  these  researches,  however,  cannot  be  considered  as  entirely  satis¬ 
factory  because  they  are  either  too  restricted  or  involve  heavy  computations .  In  his 
important  thesis  of  1942,  Nijboer  (Ref.  4)  obtained  a  simpler  and  more  satisfactory 
solution  of  the  problem  for  cases  where,  for  a  circular  aperture,  the  wave  deforma¬ 
tion  produced  by  the  geometrical  optical  aberrations  is  only  a  fraction  of  a  wavelength. 
Instead  of  using  Lommel's  series  expansions  in  terms  of  Bessel  functions,  Nijboer 
applies  expansions  of  the  aberration  function  in  terms  of  the  circular  polynomials 
introduced  by  Zernike  (Ref.  5) .  He  obtains  a  formally  different  but  equivalent  repre¬ 
sentation  to  that  of  Lommel  for  the  three-dimensional  intensity  distribution  near  the 
focus  of  an  aberration- free  optical  system.  However,  by  this  procedure,  he  is  also 
able  to  incorporate  the  effects  of  the  aberrations . 
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The  literature  on  diffraction  by  rectangular  apertures  in  the  presence  of  aberrations 
is  rather  sparse.  On  the  other  hand,  because  of  their  specific  structure,  diffraction 
patterns  of  rectangular  apertures  are  believed  more  advantageous  than  those  of  cir¬ 
cular  or  annular  apertures  in  cases  where  the  geometrical  image  location  of  an  off- 
axis  object  is  to  be  determined  by  a  two-dimensional  scanning  procedure.  For  this 
reason,  it  is  attempted  in  this  paper  to  obtain  a  new  solution  of  the  Fresnel-Kirchhoff 
diffraction  integral  for  a  stationary  rectangular  aperture .  This  solution  represents 
the  pattern  near  the  focus  of  the  optical  sys  tern  in  the  presence  of  geometrical  optical 
aberrations  due  to  a  point-source  type  object  at  infinite  distance  in  an  off-axis  position. 
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Section  2 

ASSUMPTIONS  AND  RESTRICTIONS  IN  APPLYING  THE  FRESNEL- 
KIRCHHOFF  FORMULA  TO  DIFFRACTION  PROBLEMS 


It  has  been  pointed  out  earlier  that  the  Fresnel-Kirchhoff  formula  can  be  considered 
only  as  an  approximation.  The  reason  for  this  can  be  seen  in  the  fact  that  diffraction 
theory  cannot  be  treated  in  a  completely  general  manner,  i.e. ,  it  cannot  be  rigorously 
regarded  as  describing  the  diffraction  in  telescopic  optical  systems  under  all  condi¬ 
tions.  Therefore,  it  is  deemed  advisable  to  describe,  in  this  paper,  the  assumptions 
and  restrictions  to  which  Fresnel-Kirchhoff's  diffraction  integral  is  subjected. 

Although  known  for  centuries  (Leonardo  da  Vinci  first  mentioned  the  appearance  of 
diffraction  phenomena) ,  diffraction  could  not  be  explained  on  the  basis  of  the  corpus¬ 
cular  theory  of  li^t  until  the  early  nineteenth  century.  Fresnel,  in  1816,  first  showed 
from  wave  theory  considerations  that  diffraction  phenomena  can  be  explained  by  the 
application  of  both  the  Huygens' principle  of  secondary  wavelets  and  Young's  inter¬ 
ference  principle.  Today,  the  treatment  of  diffraction  problems  is  based  on  the  con¬ 
cept  of  the  electromagnetic  wave  nature  of  light,  and  consequently  Maxwell's  equations 
are  employed  as  the  fundamental  basis  for  this  treatment . 

In  its  usual  formulation,  diffraction  theory  presumes  that  a  scalar  wave  equation  can 
be  applied  in  solving  the  problem  of  describing  the  vector  quantities  of  the  electro¬ 
magnetic  field.  From  this  point  of  view,  first  of  all,  it  has  to  be  shown  that  the  use 
of  a  sinfdc  scalar  wave  instead  of  a  vectorial  wave  is  justified  in  treating  diffraction 
problems  of  instrumental  optics.  An  example  of  this  justification  may  be  found  in 
Born- Wolf's  Principles  of  Optics ,  p.  386  (Ref.  1). 

Mathematically,  then,  the  inhomogeneous  scalar  wave  equation,  representing  a  sig¬ 
nificant  simplification  of  the  problem  to  be  solved,  is  separated  in  a  well-known 
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manner  into  two  equations  by  assuming  the  general  solution  in  the  form  of  a  product 
of  two  functions,  one  of  which  is  time  dependent  and  the  other  space  dependent.  The 
space-dependent  portion  of  the  solution,  satisfying  Helmholtz’s  homogeneous  wave 
equation,  enters  the  Helmholtz-Kirchhoff  integral  theorem  which  expresses  this  solu¬ 
tion  at  any  point  in  the  field  in  terms  of  the  values  of  the  solution  and  its  first  deriva¬ 
tive  along  the  inward  normal  at  all  points  on  an  arbitrary  closed  surface  surrounding 
the  observation  point.  It  is  to  be  noticed  that  in  order  to  obtain  the  solution,  the 
determination  of  the  appropriate  Green’s  function  is  possible  only  in  the  simplest 
cases,  i.e. ,  where  the  surrounding  surface  is  a  plane. 

Kirchhoff’s  approximation,  derived  from  the  more  general  integral  theorem,  consists 
of  a  reduction  of  the  problem  to  a  much  simpler  form. 

First,  the  closed  surface  around  the  observation  point  is  assumed  to  be  formed  by 
three  parts: 

•  A  sharp-edged  opening  in  a  plane  screen  through  which  light  from  the  object 
point  outside  the  surface  penetrates  to  the  observation  point.  If  there  is  no 
other  obstacle  between  object  point  and  diffracting  aperture,  the  wavefront 
of  the  incident  light,  of  course,  is  spherically  shaped 

•  A  portion  of  the  nonilluminated  side  of  the  screen 

•  A  portion  of  a  sphere  of  sufficiently  large  radius  with  center  at  the  observa¬ 
tion  point 

It  is  furthermore  assumed  that  the  linear  dimensions  of  the  diffracting  aperture  are 
large  compared  with  the  wavelength  of  the  monochromatic  (or  at  least  quasi- 
monochromatic)  light  emerging  from  the  source  outside  the  closed  surface  and  small 
compared  with  the  distances  of  the  object  and  observation  points  from  the  plane  screen. 

Nextt  it  is  supposed  that  the  disturbance  function  and  its  first  derivative  along  the 
inward  surface  normal  at  any  point  of  the  opening,  except  points  in  the  immediate 
neighborhood  of  the  aperture  rim,  will  not  differ  appreciably  from  the  values  of  the 
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same  quantities  encountered  in  the  absence  of  the  screen.  In  evaluating  the  diffraction 
integral,  the  contributions  from  these  quantities  encountered  at  excepted  points  are 
assumed  to  be  ne^igibly  small.  It  is  also  supposed  that  contributions  will  be  approxi¬ 
mately  zero  from  points  on  the  nonilluminated  side  of  the  screen.  In  mathematical 
formulation  these  assumptions  constitute  Kirchhoff’s  boundary  conditions.  The  addi¬ 
tional,  physically  reasonable  assumption  of  not  strictly  monochromatic  waves,  then, 
makes  the  integral  over  the  spherical  portion  of  the  closed  surface  vanish  as  the  radius 
of  the  sphere  goes  to  infinity  (see  Born- Wolfes  textbook  (Ref.  1)  fox  details).  Thus, 
the  only  remaining  portion  of  the  closed  surface  contributing  to  the  disturbance  at  the 
observation  point  is  the  aperture  itself. 

The  introduction  of  Kirchhoff^s  boundary  conditions  has  proved  satisfactory  in  many 
cases.  As  a  theory,  however,  it  does  not  indicate  the  cases  for  which  the  conditions 
are  no  longer  applicable. 

Evidently,  the  aperture  in  the  plane  screen  may  be  replaced  by  an  open  portion  of  any 
other  surface  as  long  as  the  rim  of  the  new  opening  represents  the  intersection  of  the 
new  surface  with  the  boundary  surface  of  the  pencil  of  rays  defined  by  the  object  point 
and  the  plane  aperture.  Of  course,  the  closed  surface  around  the  observation  point, 
now,  has  to  be  completed  by  the  portion  of  this  boundary  surface  between  the  edges  of 
the  plane  aperture  and  the  rim  of  the  new  opening.  The  disturbance  at  the  observation 
point  is  not  affected  if  the  contributions  to  the  integral  from  that  boundary  surface 
portion  are  negligibly  small.  This  assumption  will  certainly  hold  as  long  as  the  radius 
of  curvature  at  any  point  of  the  new  opening  is  sufficiently  large,  which  is  fulfilled  if 
one  chooses  an  arbitrary,  in  particular  a  spherical,  wavefront  as  the  new  intersecting 
surface . 

These  considerations  substantially  simplify  the  general  Fresnel-Kirchhoff  diffraction 
integral.  Its  evaluation,  then,  yields  automatically  the  discrimination  between  Fresnel 
and  Fraunhofer  diffraction.  The  simpler  case  of  Fraunhofer  diffraction,  involving 
plane  waves  emerging  from  a  point  source  at  infinity  and  diffracted  into  an  observation 
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plane  at  infinity,  is  much  more  important  in  practical  optics  than  that  of  Fresnel  dif¬ 
fraction,  which  primarily  deals  with  spherical  waves.  Unfortunately,  infinite  distances 
of  object  and  observation  points  from  the  screen  are  not  realizable  in  practice.  For 
simulating  the  ideal  case,  one  has  to  use  a  collimating  optical  system  between  the 
radiation  source  (of  finite  extension)  and  the  diffracting  aperture  in  order  to  generate 
an  approximately  plane  wave  front,  as  well  as  an  imaging  system  behind  the  screen  in 
order  to  produce  a  converging  beam  of  light  with  which  to  observe  the  diffraction  pat¬ 
tern  at  finite  distance  from  the  screen.  Nevertheless,  the  fact  that  neither  ideal  point 
sources,  even  without  using  a  collimating  system,  nor  sharp-edged  apertures  are 
realizable  in  practice,  has  to  be  considered  when  applying  the  Fresnel-Kirchhoff  dif¬ 
fraction  integral  insofar  as  boundary  effects  are  concerned. 

As  long  as  one  can  neglect  the  restrictions  and  one  stays  in  the  validity  range  of  geo¬ 
metrical  optics,  one  obtains,  in  the  Gaussian  focal  plane,  as  the  image  of  an  arbitrary 
point  source,  the  characteristic  diffraction  pattern  as  described  by  the  Fresnel- 
Kirchhoff  integral  in  the  absence  of  aberrations. 

If  one  considers  the  three-dimensional  intensity  distribution  of  the  diffraction  pattern 
near  the  focus,  one  must  extend  the  integrand  in  Fresnel-Kirchhoff ’s  integral  (as  done, 
for  example,  by  Lommel),  even  in  the  case  of  an  aberration-free  optical  system. 

As  one  exceeds  the  limits  of  Gaussian  optics,  the  wave  fronts  in  the  converging  beam 
behind  the  aperture  deviate  more  and  more  from  their  spherical  shape  because  of  the 
geometrical  optical  aberrations.  Consequently,  considerable  image  deformations  are 
to  be  observed.  If  the  aberrations  are  large  enough,  the  defect  of  the  optical  images 
obtainable,  in  general,  is  studied  by  geometrical  optics.  If,  however,  as  in  the  case 
considered  in  this  paper,  the  aberrations  are  sufficiently  small,  the  image  deformation 
has  to  be  treated  by  the  diffraction  theory  of  aberrations. 
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Section  3 

THE  DIFFRACTION  INTEGRAL  FOR  A  RECTANGULAR  APERTURE 

The  Fresnel-Kirchhoff  diffraction  integral,  in  the  specific  case  treated  in  this  paper, 
is  to  be  used  in  considering  the  disturbance  U(  P)  at  a  typical  point  P  in  the  neigh¬ 
borhood  of  the  geometrical  image  point,  P*  ,  of  a  point  P^  in  the  object  space.  It 
is  derived  from  the  geometrical  configuration  encountered  in  Fig.  1. 

The  planes  that  are  perpendicular  to  the  optical  axis  of  the  image-forming  instrument 
and  contain  the  points  P^  or  P*  are  usually  denoted  as  object  or  image  planes, 
respectively.  The  lateral  magnification  factor  M  between  these  planes  is  determined 
by  the  classical  rules  of  geometrical  optics.  The  plane  rectangular  aperture  acts  as 
the  entrance  pupil  of  the  optical  system.  Application  of  the  same  classical  rules  yields 
the  exit  pupil  location  and  the  Gaussian  lateral  magnification  factor  M’  between  both 
pupil  planes  which  are  also  considered  to  be  perpendicular  to  the  optical  axis.  The 
meridional  plane,  defined  by  P^  (and  P*)  and  the  optical  axis,  is  taken  as  a  refer¬ 
ence  plane. 

Two  Cartesian  coordinate  systems  are  introduced  with  origins  at  C  and  P*  ,  respec¬ 
tively.  The  positive  Z-axes  have  the  direction  of  the  ray  from  C  to  P*  (direction 
of  light  propagation);  the  other  axes  are  mutually  parallel  and  have  equal  positive 
directions . 

For  convenience,  the  Y-axes  are  chosen  in  the  above-defined  meridional  plane.  Points 
in  the  system  with  origin  at  C  are  described  by  coordinates  x  ,  y  ,  z  ,  and  points  in 
the  other  system  by  X  ,  Y  ,  Z  . 

The  aperture  (in  the  exit  pupil  plane)  is  now  replaced  by  that  portion  of  the  Gaussian 
reference  sphere  centered  at  which  passes  through  C  and  approximately  fills 
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this  aperture.  The  actual  wave  front  passing  through  C  ,  in  general,  will  deviate 
slightly  from  the  spherical  shape  because  of  the  small  aberrations  being  assumed.  An 
arbitrary  ray  diffracted  into  the  image  space  in  the  direction  of  the  observation  point 
P(X  ,  Y  ,  Z )  in  the  neighborhood  of  P*  passes  the  aperture  at  Pj(X^  .  )  and 

intersects  the  Gaussian  reference  sphere  and  the  actual  wave  front  at  Q(X’  ,  Y*  ,  Z’) 
and  Q(X’  ,  Y’  ,  Z^ ) .  It  also  is  possible  to  describe  Q  by  its  coordinates  x  ,  y  ,  z  . 

Of  course,  these  coordinates  may  be  deduced  from  the  other  ones  by  the  transformation 

X^  -  x;  Y^  =  y  ;  Z'  =  z  -  R 

Because  P  is  a  point  in  the  immediate  neighborhood  of  P*  and  the  aperture's  linear 
dimensions  are  assumed  to  be  small  compared  to  the  radius  R  of  the  Gaussian 
reference  sphere,  the  angle  x  between  the  diffracted  ray  and  the  normal  to  the  sphere 
at  Q  will  be  small.  Thus,  it  can  be  assumed  that  the  variation  of  the  inclination 
factor  in  Fresnel -Kirchhoff's  general  diffraction  integral  will  become  small  as  Q 
explores  the  portion  of  the  reference  sphere  corresponding  to  the  aperture  area.  For 
this  reason,  this  variation  may  be  neglected. 

Then,  denoting  the  distance  between  a  typical  point  Q  on  the  reference  sphere  and  the 
observation  point,  P,  by  s  ,  and  remembering  that  the  distance  between  Q  and  Q 
defines  the  wave  aberration  ^  ,  the  diffraction  integral  takes  the  form 

A  f(  exp  [jk(4>  +  s  -  R)] 

U(P)  =  -J  t  JJ  ■  — L  CIS 


where: 

k  =  27r/X 

X 

A 

dS 


the  wave  number 

the  wavelength  of  the  monochromatic  or  quasi -monochromatic  light 

the  amplitude  of  the  wave  (assumed  to  be  substantially  constant 
over  the  wave  front) 

the  surface  element  on  that  portion  of  the  reference  sphere  over 
which  the  integration  is  to  be  extended 
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To  obtain  an  approximate  solution  of  the  integral,  for  observations  points  P  in  the 
region  of  the  geometrical  image  point  ,  the  quantity  s  in  the  denominator  of  the 
integrand  usually  is  replaced  by  the  radius  R  of  the  reference  sphere .  The  surface 
element  is  taken  as 


dS  =  dx  dy 

Taking  the  integration  limits  as  ±M^a  and  ±M’b  ,  where  2a  and  2b  are  the  linear 
dimensions  of  the  entrance  pupil  (diffracting  aperture),  is  restricted  to  cases  where 
only  relatively  small  incidence  angles  of  the  principal  ray  are  involved  (supposition 
of  small  aberrations). 


From  the  analytical  expressions  for  the  distance  s  ,  the  radius  R  of  the  reference 
sphere  and  the  distance  r  between  P*  and  P  ,  supposed  to  be  small ,  one  obtains 


s  -  R 


In  this  expression,  all  small  terms  of  second  order  have  been  neglected. 
Hence 


M^a  M*b 


U(X,  Y,  Z)  -  -j 


A  jkZ 
\  ® 


m  a  AVI  u  t 

J  J  exp  jk  - 

M»a  -M'b  L  \ 


Xx  + 


R 


15^  ^  Z  xl 
2 


R 


■) 


— 9  dxdy 
R 


For  convenience,  new  integration  variables  are  introduced,  and  new  parameters 
describing  the  observation  point  location  in  the  geometrical  image  region  are  defined. 
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Integration  variables: 

X  =  IVt  = 

dx  =  R  d| 

Integration  limits; 

X  =  ±M'a 

=  ±M'a/R  =  ±Mq 

4  =  ±1 

Optical  coordinates; 


y  =  Rj^  =  Rv^tj 

dy  =  R  d^'  =  Rt*  dn 
o 


y  =  ±M'b 
V  =  ±M'b/R  = 

TJ  =  ±  1 


u  =  kZ 


V  =  kX 


w  =  kY 


Hence 


1  1 

U(u,v,w)  =  if  Jexph(k*-|Mo4^-VMo^  dTJ 

-1-1  *■ 

In  the  case  of  an  aberration- free  optics,  i.e,,  for  ^  0  ,  obviously  the  two- 

dimensional  diffraction  integral  is  expressible  as  the  product  of  two  one-dimensional 
integrals .  It  seems  appropriate  at  this  point  to  attempt  an  adequate  separation  also 
for  the  optical  system  affected  by  aberrations. 

Usually,  to  obtain  an  analytical  expression  for  the  wave  aberration,  two  mutually 
parallel  and  equally  oriented  Cartesian  coordinate  systems  are  introduced  with 
origins  on  the  optical  axis  of  the  image- forming  instrument  at  the  centers  of  the 
object  and  image  planes,  respectively.  The  positive  Z-axes  are  taken  along  the 
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optical  axis  in  direction  from  the  object  to  the  image  space.  The  object  and  image 
points  are  defined  by  P^(X^,y^,0)  and  P|(Xj,Yj,0),  respectively.  If  M  is 
the  Gaussian  lateral  magnification  factor  between  the  object  and  image  planes,  one  has 

X?  =  MX  and  Y*  =  MY 

1  o  1  o 

(Arbitrarily  placing  the  Y-axes  into  the  meridional  reference  plane,  defined  in  the 
beginning,  one  will  obtain  =  X|  =  0 .) 

The  center  of  the  exit  pupil  plane  is  denoted  by  C(0  ,  0  ,  -  Dj)  where  Dj  is  the  dis¬ 
tance  between  the  exit  pupil  and  image  planes.  Again,  the  Gaussian  reference  sphere 
is  centered  at  P ^  and  passed  through  C  .  The  actual  wave  front  also  passes  through 
C  •  If,  for  the  moment,  diffraction  effects  are  omitted,  a  specific  nondiffracted 
ray  emerging  from  will  meet  the  entrance  pupil  at  *  ^o^  ’ 

the  exit  pupil  at  ,  Z^) ,  and  the  image  plane  at  ,  Y^^ ,  0) , 

(see  Fig,  2) . 

If  the  ray  [  P^  P^  ]  intersects  the  actual  wave  front  and  the  Gaussian  reference 
sphere  at  Q^{X' ,  Y*,  Z')  and  Qj^(X’,  Y^  Z’)  ,  respectively,  the  aberration  func¬ 
tion  is  defined  as  the  distance  between  these  two  points 

^  =  [Q^Qll 

By  making  use  of  the  fact  that  both  and  C  are  located  on  the  same  wave  front 
one  obtains 

*  =  -  IPoCl 

i.e.,  $  is  expressible  as  the  difference  of  two  light  paths.  Because  Q  lies  on  the 
Gaussian  reference  sphere ,  its  coordinate  Z  can  be  eliminated  if  both  optical  paths 
are  expret.3ed  in  the  notations  of  the  point  characteristic,  and  one  obtains  finally 

$  =  *(Xg,Y^;  X',Y)  (Ref.  1,  p.  202) 
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This  function  is  known  to  be  expressible  in  expanded  form  by 

$  =  $(4)  +  *(6)  +  *(8)  ^ 

where  is  a  polynomial  of  degree  2k  in  the  four  coordinates  X^,Y^,X',Y', 

containing  these  quantities  only  in  combinations  of  powers  of  the  three  scalar  invariants 

(X^  +  Y^)  .  (X'^  +  Y'2)  ,  (X^X'+  Y^Y') 

If  diffraction  is  involved,  the  specific  ray  taken  into  consideration  is  diffracted  at  the 
above  defined  point  >  2^)  directions  of  the  space  behind  the  screen. 

If  the  diffraction  effect  is  to  be  observed  at  point  P(Xj^ ,  in  the  neighborhood 

of  P*,  the  corresponding  diffracted  ray  emerging  from  P^  will  intersect  the  exit 
pupil  at  P’  (X^  ,  )  .  The  intersection  points  of  the  ray  [  P’  P  ]  with  the  actual 

wave  front  and  the  Gaussian  reference  sphere  may  now  be  defined  by  Q(X*,  Y’,  Z*) 
and  Q(X*,  Y*,  Z* ) ,  respectively.  The  wave  aberration  again  is  given  as  the  distance 
between  these  points: 

^  =  [QQ]  =  [PqQ]  -  [PqC]  =  ^(Xo,Y^;  X\Y') 

The  general  expression  obtained  for  the  aberration  function,  therefore,  remains 
unchanged. 

According  to  the  assumption  of  the  smallness  of  the  ratios  between  the  exit  pupiPs 
linear  dimensions  and  the  radius  of  the  Gaussian  reference  sphere,  the  coordinates 
X*,  Y^  may  be  replaced  by  the  coordinates  X^  ,  Y^  of  the  intersection  point  P’ . 

Instead  of  the  real  coordinates,  SeidePs  variables  are  introduced  by  the  transformation 


V  -  o  Q 

n  X 
o  o 


Y 

"o  "  11  X 
o  o 


3-8 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


6-90-62-116 


X, 

DiXi 

DiYi 

1 

1 

X’ 

=  X  i 

Y'  =  X  77 

o 

o  o 

o 

o 

o' 

x'l 

= 

where 


D 

D 


n 

n 


o 

1 


the  distance  between  the  object  and  entrance  pupil  planes 

the  distance  between  the  entrance  pupil  and  image  planes 

refractive  index  of  the  object  space 

refractive  index  of  the  image  space 

length  unit  in  the  entrance  pupil  plane 

length  unit  in  the  exit  pupil  plane 


Because,  within  the  accuracy  of  Gaussian  optics,  one  has 


yi  =  ^0 


^1  =  ^0 


’’i  =  \ 


one  obtains  the  relation,  with  respect  to  the  object  and  geometrical  image  points, 


X 


o 


y 


★ 

1 


For  simplification,  may  be  replaced  by  R,  the  refractive  indices  n^  and  n^^ 
and  the  length  unit  may  be  taken  equal  to  unity.  Then,  characterizes  the 
length  unit  in  the  exit  pupil,  as  well  as  the  magnification,  M’  =  (dimensionless 

in  this  case ,  of  course)  between  the  two  pupil  planes . 

The  aberration  function  expressed  in  Seidel’s  variables  takes  the  form 

k  =  2 
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where  the  polynomial  is  a  linear  power  combination  of  degree  (2k)  in  the 

three  basic  quantities 


2  *2  ^  *2  2  .2^2  2  ^  * 

>*0  =  X*  +  y*  ;  p  =  +  J?!  :  X  ^  ^1^1 


With  respect  to  the  primary  or  third-order  aberration  one  has  in  usual  notation 


$ 


(4) 


-  I  Bp^  -  -  I  Dr2p2  +  Er^K^  +  Fp^/c^ 


The  tunction  corresponding  to  the  fifth-order  aberrations  contains  the  quantities 

6  6  24  24  42  42  42  24  ^(8)  8  8  2  6 

p  ,  K  ,  r^p  ,  r^K  ,  r^p  ,  t^k  ,  p  /c  ,  p  k  ;  'the  quantities  p  ,k  ,  r^p  , 

26  44  ^4  4  62  62  62  44  26  , 

r^K  ,  r^p  ,  T^K  ,  r^p  ,  t^k  ,  p  k  ,  p  k  ,  p  k  ,  and  so  on  to  infinity. 

It  is  now  recommended  that  Seidel  *s  variables  be  related  to  the  integpration  variables 
introduced  above.  This  yields 


xj  =  -X^x*/R 

xA^  =  Rp^lAi 


yj  -  -x^Y*/R 


’ll  y/^i  =  R‘'o^Ai 


By  this  transformation  Seidel’s  third-order  aberration  function  is  expressible  in  the 
integration  variables.  The  wave  aberration  may  be  simplified  by  arbitrarily  choosing 
=  0  .  Thus,  one  obtains  in  explicit  form 


y*^  =  (■^)  ^  constant 

.2  ^  2  /  R  \^  /  2.2  ^  2  2\ 

i"^i  -[tJ  "  v) 


y>i  = 
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and 


(4)  _ 


cY;2.y 


- 1  Bvf  .  .y) 

-  B(^j  Iff  V 
-(^f  n(.^o^v  .y) 


The  expressions  for  the  aberration  functions  of  higher  order  are  obtained  in  the  same 
way. 

It  is  seen  that  the  wave  aberration  consists  of  a  linear  combination  of  powers  in  4 
and  7} .  For  this  reason,  the  exponential  term  in  Fresnel-Kirchhoff’s  formula  con¬ 
taining  ^  in  its  exponent  cannot  be  separated  immediately  into  two  terms  containing 
only  powers  of  4  or  i? . 
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Section  4 

EXPANSION  OF  THE  INTEGRAND  IN  FRESNEL-KIRCHHOFF’S  FORMULA 

As  mentioned  in  earlier  sections  of  this  report,  the  application  of  Zernike^s  circular 
polynomials  enabled  Nijboer  to  incorporate,  in  the  case  of  diffraction  by  a  circular 
aperture,  the  effect  of  aberrations  by  expanding  the  integrand  in  Fresnel-Kirchhoff’s 
formula.  This  yielded  an  approximate,  but  nevertheless  quite  satisfactory,  solution. 
Zemike^s  circle  polynomials  are  complex  functions  in  two  real  variables  of  the  form 

(p  sin  0  ,  p  cos  e  )  =  (P)  e^”^^ 

where  m  |  0  and  n  =  0  are  integers,  n  ^  |  m  j  ,  and  n  -  i  m  j  is  even. 


They  satisfy  the  orthogonality  relation 
27rl 

/ / R“(P)R“'(p)  exp  [j(m  -  m')0]  p  dpde  = 
o  o 

where  5^  is  the  well-known  Kronecker  symbol.  This  property  is  essential  in  Nijboer’s 
integrand  expansion  because  this  expansion  necessarily  involves  the  occurrence  of 
Zernike's  radial  polynomials  rJJ^  (p)  in  all  portions  of  the  integral. 

The  rectangular  aperture,  unfortunately,  does  not  have  circular  symmetry;  rather,  it 
is  characterized  by  finite  limits.  Although  these  limits  may  be  supposed  to  be  symmet¬ 
rical  to  the  aperture  center,  they  are  different  for  each  of  the  two  variables,  x  and  y. 
The  only  adequate  polynomials  having  orthogonality  properties  in  a  finite  integration 
interval  are  seen  in  Legendre's  polynomials  (this  is  also  stated  in  a  small  note  by 
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Zemike  and  Nijboer  (Ref.  6)  -  *'La  theorie  des  images  optiques/'  p.  228).  To  become 
orthogonal  between  finite  limits  (for  example  or  respectively)  the  com¬ 

mon  expressions  for  these  polynomials  can  easily  be  modified.  However,  it  is  be¬ 
lieved  more  convenient  to  transform  the  variables  x  and  y  into  an  appropriate  form 
such  that  the  ordinary  Legendre  polynomials,  with  orthogonality  property  in  the  inte¬ 
gration  interval  ±  1 ,  can  be  used  in  expanding  the  terms  of  the  diffraction  integral. 
For  this  reason,  in  the  previous  section,  the  dimensionless  integration  variables, 

I  and  Tj ,  have  been  introduced. 

The  first  step  in  the  expansion  to  be  performed  refers  to  the  exponential  terms  con¬ 
taining  the  first  and  second  powers  of  |  and  r]  separately.  The  second  step  concerns 
the  exponential  term  containing  the  wave  aberration  ^  . 

According  to  Bauer  (Ref.  7)  one  has 


exp  ( jp  cos  6  ) 


(j)®(2s  +  1) 


s=o 


pl/2 


Pg  (cos  0) 


Substituting 


cos  $ 


2 


and 


cos  0  =  -  ^ 


and  noting  that 


P(-cos0)  =  (-l)®P(cos9) 

S  VJ 
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one  obtains 


(->  i  * 


/:  ^  (-I)N2s  *  1)  p  ,j2, 


exp  (-jvp^o 


=  ^71  i  <-!)“' 


(vPo)^ 


and  similar  expressions  for  the  exponential  terms  containing  the  variable  7]  in  their 
exponents . 

The  combined  exponential  term  exp|- j  V2  +  vp^«)]  is  expressible  in  the  form 

2 

of  a  double  sum  containing  products  P  )  P_t(| ) .  The  polynomials  involved, 

s  s 

unfortunately,  are  not  orthogonal  as  can  easily  be  shown. 


However,  the  Legendre  polynomials  may  be  defined  by 


Pod)  - 1 


where  6^  is  a  common  denominator  for  all  coefficients.  Similarly,  the  powers  of 
i  may  be  expressed  by 


-  2  ^  2  [(“  It  -  «  *  ‘)]  V2.«) 


6^  being  again  a  common  denominator  for  all  coefficients . 
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From  these  definitions,  it  is  obvious  that  by  replacing,  in  P  (|) ,  the  variable  | 

2  2  °o 

by  4  ,  i.e, ,  by  expressing  P  (|  )  explicitly  in  terms  of  |^  ,  and  by  substituting, 

s 

for  the  powers  of  |  in  the  linear  combination  obtained,  their  expansions  in  terms  of 

2 

common  Legendre  polynomials  in  f,  ,  P  (|  )  is  represented  by  a  finite  double  sum 

s 

in  terms  of  those  ordinary  Legendre  polynomials.  Moreover,  it  is  seen  that  the 
exponential  term  exp  [  -  j  72  +  is  expressible  in  the  form  of  a  quad¬ 

ruple  sum  in  terms  of  products  of  common  Legendre  polynomials  where  every  single 
term,  of  course,  has  the  usual  orthogonality  property. 

A 

The  exponential  term  e*’  ,  accordingly,  is  expanded  into  its  power  series  in  (jk  ) 


1  +  j 


JL  *  k_ 

1  *  ^  ”  a  I 


* 


k: 

3! 


+  k  *4 
+  TT  ^  ■ 


(jk«)^ 


£  =  o 


In  this  expression,  the  wave  aberration  ^  ,  in  general,  can  be  considered  as  a 
double  series  in  terms  of  power  products  where  n  =  0,l,2,...oo, 

m  =0,l,2,...n. 


00  n 

y  y  K 

Z-f  n-m ,  m 

n  =  o  m  =  0 


where  the  constants  K  _  are  dependent  on 
n-m ,  m 

•  The  off-axis  position  of  the  geometrical  image  point,  X*  and  Y* 

•  The  linear  dimension  parameters ,  and  ,  of  the  aperture 

•  The  radius  R  of  the  Gaussian  reference  sphere 

•  The  unit  length  of  the  exit  pupil  plane 

•  Certain  aberration  coefficients 
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r 

{ 

? 


{ 


Taking  the  power  of  the  aberration  function  does  not  change  the  general  expres¬ 
sion,  except  that  the  coefficients  of  the  expansion  now  are  different.  This  may  be 
expressed  by 


i  i  4  m  m 

n  =  o  m  =  o 


where 


and  yields 


JO-  00  n 

I  2<- 

i  =  o  n  =  om  =  o 


.n-m  m 


Evidently,  in  every  single  product  of  this  expression,  the  powers  of  ^  and  7j  can 
be  expanded  in  terms  of  Legendre ^s  polynomials. 

By  these  expansions  of  the  integrand  terms ,  one  can  divide  the  two-dimensional  dif¬ 
fraction  integral  into  two  one-dimensional  integrals 


if  I  expn(k.t  -  ^  -  vnj  -  I  d$  dv 

■»  *1  • 


-1-1 


oo  n 


=  1  V.  «'■>'  2  2 


f  =  o 


n  =  o  m=o 
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Then,  the  disturbance  at  the  observation  point  P  is  given  by 


00  n 


Jf=0 


n  =  o  m=o 


where 


+  n  2s'  +  1  f 


B^0  8’-0 


(VPo)*'“  -1 


^  \  ^  I  n  '^t+l/2(“‘'oA)  2t'+l  '^t'+l/2^'^V  f  _m„  .  2,„  _  ^ 


t-O  t*-0 


(^v  -1 
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Section  5 

EVALUATION  OF  THE  DIFFRACTION  INTEGRAL 


After  the  reduction  of  Fresnel-Kirchhoff’s  diffraction  formula  for  the  rectangular 
aperture  to  a  much  simpler  form,  the  remaining  integrations  to  be  performed  con¬ 
cern  only  integrals  of  the  type 


/ 

- 1 


x*^Pi.(x^)Pp,(x)  dx 


where  x  stays  for  any  of  the  variables  ^  and  rj;  r  ,  r’  replace  any  of  the  subscripts 
s  ,  s*  and  t ,  t* ,  respectively;  and  k  replaces  any  of  the  exponents  (n-m)  and  m. 
In  order  to  evaluate  this  integral ,  obviously  one  will  have  to  make  use  of  the  orthogo¬ 
nality  property  of  Legendre^s  polynomials 


/ 
- 1 


P^(x)P^(x)dx  = 


6^  being  the  Kronecker  symbol , 


It  is  deemed  advisable  to  introduce  new  polynomials  defined  by 


izm 

2 


iim 
2 


p  =  o 


p  =  o 


for  r  ?  1 
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Expres$ing  the  powers  of  x  by  their  expansions  in  terms  of  Legendre  polynomials, 


2r-Ht-4p 


°‘2r+k-4p,K 


^2r+k-4p  ^2r4-k-4p-2x<’‘) 


2r+k-4p-(l) 


o2r+k-4p  ^  rt-2#c  4r  +  2k  -  8p  -  4k  +  1 
^  ^  2r  +  k  -  4p  -  2k  +  1 

K  =0 

^  [(2r.k-4p-«)/(4r.^-  8.-^2,.  1)], 


^2r+k-4p-2x^*) 


k-(l) 

2 


6^  ^k-2K^^^ 


one  obtains 


r-(l)  2r+k-4p("l) 

,  y  ^  y  °2r-Hc-4p,« 


/’4) 


V  “k.K 

Z  -6-  V2k(") 


The  values  of  ,  /3.  ,  6.  ,  6!  are  known  from  Salzer's  (Ref.  8)  tables  of  the  powers  of 
X  in  terms  of  Legendre  polynomials  (up  to  order  twenty-four)  and  from  Tallquist's 
(Ref.  9)  tables  of  the  Legendre  polynomials  (up  to  order  sixteen)  and  can  easily  be 
extended  to  higher  orders. 
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The  expression  obtained  for  the  polynomials  /^(x)  may  be  more  simplified  by  taking 
into  account  the  fact  that  some  of  the  coefficients  of  the  expansion  are  to  be  related  to 
a  Legendre  polynomial  of  the  same  order.  This  yields 

2r4k-(l) 

^r^*^  "  2  ^2r+k-2l(*) 

i  =  o  ^r 

where  D?'  is  a  common  denominator  for  all  the  coefficients.  The  coefficients  j 
and  have  been  computed  for  OS  r  S  12  and  0  S  k  S  24  with  the  restriction 
[2r  +  k-  (l)]/2S  12,  and  are  available  by  request  (Ref.  10). 


By  this  procedure,  the  evaluation  of  the  typical  integral  is  reduced  to 


X  .4. 

j  =  J/^r(x)Pj.,(x)  dx 


2r4k-(l) 


2 


i  =  o  ®r  “  ^ 


2r4k-(l) 

- 2- - ..k 


V  Nr.i  2 
Z  ok  2r'  +  1 


i  =  o  r 


where  necessarily  r'  =  2r  +  k  -  2i . 
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For  convenience  the  functions  are  defined 


/  2r'  +  1  '^r'+l/2^^^ 

=  <"J>  2  ITT 


2r+k-(l) 
2 


-1  i  =  o 


iC 

■-i  ^r,i  •^2r-ric-2i+l/2^^^ 

„k  ttt 


where  z  stays  for  any  of  the  variables  ( )  and  ( ) . 


Substituting  these  results  into  the  expressions  for  and  W™  obtained  above,  one 

has  finally 


S  =  0 


1) 


Wj“(u,  w) 


t  =  o 


1) 


where 


2s+n-m-(l) 
2 


G"-“(vpo) 


2 


a=o 


J2s+n-m-2a+l/2(vfio) 

(- 1)  ^ — "172 - 


D 


n-m 

s 


(vPq) 


2t+m-(l) 
2  ■ 


g“(w.o)  = 


t-r  Nt?T  J2t+in-2T+l/2('^^) 

^ — ZTjrt — 


T  =  0 


(wi^q) 
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Normally ,  the  three-dimensional  intensity  distribution  is  now  obtained  as 


I(u,v,w)  =  |U(u,v,w)r  =  I, 


00  n 


£  =  o 


n  =  o  m  =  o 


where 


I  =  (4  V 

o  \X  Of 


represents  the  intensity  at  the  center  of  the  aberration- free  diffraction  pattern. 


The  solution  in  the  Gaussian  image  plane  is  of  special  interest.  Because  of 


lim 


1/2 


=  /I  6° 

V  TT  r 


one  obtains 


n-m-(l) 
2 


Vj  (O.v) 


wJ"(o,w) 


.  .,n-m /iF  ^n-m,  v  ^  .vU-m  /tF  ^  ,,o-^o,cr 

H)  (v^o)  =  H)  V2  z 


(T^.cr  ’^n-m-2o’+l/2^'^^o^ 

172 


0  =  0 


(VMo)^ 


m-(l) 
2 


.  ..m  fit  \  .  -.T  ^o,  T  ‘^m-2T+l/2^^^o^ 

'  <-'»  h  2  '-*> 

T  =  O  O  'o' 
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Hence 


1(0,  v,w)  =  I 


00  n 


I  n  <»'  I  2  I 


!  =  o  n=o  m  =  o 


n-m-(l)  m-(l) 

V  ^n-m-2CT^  ^  ^o,  t '^m-2T+l/2^'''*^o^ 


0  =  0 


D 


n-m 


(vMq) 


V2 


V  t”o 
D 


7  =  0 


m  ~  .1/2 

1  (wi^  ) 

o  'o' 


This  yields  for  the  aberration-free  optical  system,  characterized  by  n  =  m  =  0  , 


1(0,  v,w)  =  I 


TT  Jl/2<^^)  Jl/2<^V 


/sin  v/x  X^^sin 

)  \  ^"o  / 


i.e. ,  the  classical  expression  for  the  diffraction  pattern  due  to  a  rectangular  aperture. 


Furthermore,  the  three-dimensional  diffraction  pattern  intensity  distribution  near  the 
focus  of  an  aberration-free  optical  system  is  also  believed  to  be  of  interest.  This 
distribution  is  represented  by: 

2  2 

I(u.v.w)  =  Io|v°(u.v)|  |w°(u,w)| 


where 


V°(u.v) 


f  ^  (-j)®(2s  +  1) 
s  =o 


Gs(v.Mo) 
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W°(u.w) 


i  i  »■-  •  •  ^3®  •i-v 


t  =  0 


(u^2/2)1 


G^(vu  ) 
s' 


D 

=  2 


(-ly 


5-CT  Ng  J2s-2g+l/2<VMo) 


a  =  o 


(v^o) 


1/2 


G°(w.o) 


^  /  I.t-T  ^t,T  ^2t-2T^-l/2<^''o) 

^  D° 

T  =  O  t  O 


This  equation,  of  course,  yields  the  above-derived,  classical  expression  for  the  dif¬ 
fraction  pattern  of  a  rectangular  aperture  when  taking  u  =  0,i.e.,s=t=0. 
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Section  6 
SUMMARY 


The  Fresnel-Kirchhoff  integral  at  an  arbitrary  observation  point  in  the  neighborhood  of 
the  geometrical  image  of  a  specific  object  point  represents  the  disturbance  due  to  the 
diffraction  by  a  rectangular  aperture  in  the  presence  of  geometrical  aberrations.  In 
this  paper  it  is  shown  that  this  integral  can  be  expanded  into  a  series.  The  expansion 
is  based  primarily  on  the  application  of  Bauer’s  formula  containing  Bessel  functions  of 
the  order  (n  +  1/2)  where  n  is  an  integer  and  also  on  the  representation  of  the  expo¬ 
nential  term  containing  the  wave  aberration  in  its  exponent  by  a  power  series.  By  this 
procedure ,  a  significant  simplification  is  obtained  which  is  seen  in  the  conversion  of  the 
two-dimensional  diffraction  integral  into  a  summation.  The  single  terms  are  divided 
into  products  of  two  functions  where  only  one -dimensional  Integrations  are  to  be 
performed. 

To  obtain  a  general  solution  of  the  diffraction  integral,  new  polynomials  are  introduced 
and  expanded  in  terms  of  Legendre  polynomials.  Applying  the  orthogonality  property 
of  Legendre’s  polynomials  yields  the  definition  of  special  functions,  containing  essen¬ 
tially  linear  combinations  of  Bessel  functions  of  the  order  (n  +  1/2). 

By  these  special  functions,  the  disturbance  and  the  three-dimensional  intensity  distri¬ 
bution  about  the  arbitrary  observation  point  in  the  neighborhood  of  the  geometrical 
image  point  are  expressible  in  a  general  form. 

Particular  solutions  are  obtained  for  the  two-dimensional  intensity  distribution  of  the 
diffraction  pattern  in  the  Gaussian  image  plane  in  the  presence  of  aberrations  and  for 
the  three-dimensional  intensity  distribution  near  the  focus  of  an  aberration-free  optical 
system. 


6-1 


LOCKHEED  MISSILES  de  SPACE  COMPANY 


6-90-62-116 


Section  7 
REFERENCES 


CITED  REFERENCES 

1.  M.  Born  and  E.  Wolf,  Principles  of  Optics ,  'Pergamon  Press ,  New  York,  1959, 
p.  374,  alsopp.  76  and  202. 

2.  E.  C.  J.  V.  Lommel,  Abh.  d.  Bayerischen  Akademie  d.  Wissenschaften ,  Vol.  15, 
Munchen,  Abt.  2,  1885,  p.  233;  Abt.  3,  1886,  p.  531 

3.  E.  Wolf,  Reports  on  Progress  in  Optics ,  Vol.  14,  Physical  Society,  London, 

1951,  p.  95 

4.  B.  R.  A.  Nijboer,  Thesis,  University  of  Groningen,  1942 
Physica,  Vol.  10,  1943,  p.  679;  Vol,  13,  1947,  p.  605, 

5.  F.  Zemike,  Physica ,  Vol.  1,  1934,  p.  689 

6.  F.  Zernike  and  B.  R.  A.  Nijboer,  "La  thdorie  des  images  optiques,"  Editions 
de  la  Revue  d'Optique,  Paris,  1949,  p.  227 

7.  Bauer,  J.  fiir  Mathematik,  Vol.  56,  1859,  p.  104 

8.  H.  E,  Salzer,  Math.  Tables  and  Other  Aids  for  Computation,  Vol.  3,  1948/49, 

p.  16 

9.  Hj.  Tallquist,  Act.  Soc.  Fennicae,  Nova  Series  A,  Vol.  2,  No.  4,  1937,  p.  1 

10.  H.  P.  Greinel,  Functions  and  Tables  of  Coefficients  Needed  for  the  Description 
of  the  Diffraction  Pattern,  in  the  Presence  of  Aberrations,  Due  to  a  Stationary 
Rectangular  Aperture,  Internal  Report,  Electronic  Sciences  Lab.,  Lockheed 
Missiles  &  Space  Compai^,  Palo  Alto,  Calif. 


7-1 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


6-90-62-116 


NOT  SPECIFICALLY  CITED  REFERENCES 

1.1.  H.  V.  Helmholtz,  J.  ftir  Mathematik.  Vol.  57,  1859,  p.  7 

1.2.  G.  Kirchhoff,  Berliner  Berichte,  1882,  p.  641;  Annalen  der  Physik.  (2), 

Vol.  18,  1883,  p,  663 

1.3.  H.  A.  Buchdahl,  Optical  Aberration  Coefficients ,  Oxford  Univ.  Press,  1954 

1.4.  M.  Herzberger,  Modern  Geometrical  Optics ,  New  York,  Interscience,  1958 

1.5.  L.  Seidel,  Astronomische  Nachrichten .  Vol,  43,  No.  1027,  1856,  p.  289; 

No.  1028,  p.  305;  No,  1029,  p.  321 

1.6.  K.  Schwarzschild,  Abh.  d.  Koeniglichen  Gesellschaft  d.  Wissenschaften. 
Goettingen,  Math.-phys.  Klasse,  Vol.  4,  No.  1,  2,  and  3,  1905/1906 

1.7.  A,  Kohlschuetter,  Dissertation,  Goettingen,  1908 

1.8.  G,  C.  Steward,  Trans,  of  the  Cambridge  Phil.  Soc. ,  Vol.  23,  1926,  p.  235 

1.9.  M.  Herzberger,  J.  Optical  Soc.  Am.,  Vol.  29,  1939,  p.  395 

1.10.  F.  Wachendorf,  Optik,  Vol.  5,  1949,  p.  80 

1.11.  J.  Focke,  Dissertation,  Leipzig,  1951 

1.12.  N.  Chako,  Trans.  Chalmers  Univ.  of  Technology,  Gothenburg,  1957,  p.  191 

1.13.  H.  A.  Buchdahl,  J.  Optical  Soc.  Am. .  Vol.  46,  1956,  p.  941;  Vol.  48,  1958, 
pp.  563,  747,  757;  Vol.  49,  1959,  p.  1113;  Vol.  50,  1960,  pp.  534,  539,  678 
Vol.  51,  1961,  p.  608 

4.1.  K.  Nienhuis,  Thesis,  University  of  Groningen,  1948 

4.2.  K.  Nienhuis  and  B.  R.  A.  Nijboer,  Physica,  Vol.  14,  1948,  p.  590 

5.1.  F.  Zemike  and  H.  C.  Brinkman,  Proc.  of  Section  of  Sciences,  Akademie  van 
Wetenschappen,  Amsterdam,  Vol.  38,  1935,  p.  161 

5.2.  A.  B.  Bhatia  and  E.  Wolf,  Proc.  of  the  Cambridge  Phil.  Soc. .  Vol.  50,  1954,  p.  40 

8.1.  G.  N.  Watson,  A  Treatise  on  the  Theory  of  Bessel  Functions.  Cambridge  Univ. 
Press,  1958,  p.  128 


7-2 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


